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ABSTRACT 

We present the effects of the vertical convection on the structure and luminos¬ 
ity of the neutrino-dominated accretion flow (NDAF) around a stellar-mass black 
hole in spherical coordinates. We found that the convective energy transfer can 
suppress the radial advection in the NDAF, and that the density, temperature 
and opening angle are slightly changed. As a result, the neutrino luminosity and 
annihilation luminosity are increased, which is conducive to achieve the energy 
requirement of gamma-ray bursts. 

Subject headings: accretion, accretion disks - black hole physics - convection - 
gamma-ray burst: general - neutrinos 


1. Introduction 

The central engine of gamma-ray bursts (GRBs) is usually modelled as the system 
consisting of a rapidly spinning stellar black hole and an extremely optically thick disk 
with high density and temperature, namely neutrino-dominated accretion flow (NDAF). The 
neutrino radiation is the main cooling mechanism, which possibly has the ability to power 
GRBs via neutrino annihilation outside the disk. A lot of works on this model had been done 
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in the past decades, which can be divided into the following categories: the researches on the 
stable or time-dependent radial characteristics of the disk, such as the dynamics, structure, 
components, and luminosity of the NDAF model (see, e.g., Popham et al. 1999; Narayan 
et al. 2001; Di Matteo et al. 2002; Kohri & Mineshige 2002; Kohri et ah 2005; Janiuk et 
ah 2004, 2007; Lee et ah 2005; Gu et ah 2006; Chen & Beloborodov 2007; Kawanaka & 
Mineshige 2007; Liu et ah 2007; Kawabata et ah 2008; Zhang & Dai 2009; Kawanaka et ah 
2013; Li & Liu 2013; Xue et ah 2013) or magnetic NDAF model (e.g., Lei et ah 2009; Janiuk 
& Yuan 2010; Luo et ah 2013), and the corresponding applications to GRBs observations 
(e.g., Reynoso et ah 2006; Lei et ah 2007; Lazzati et ah 2008; Liu et ah 2010b, 2012b; 
Romero et ah 2010; Kawanaka & Kohri 2012; Sun et ah 2012; Hou et ah 2014a, b), on the 
vertical characteristics of the model, and the applications to observations (e.g., Sawyer 2003; 
Liu et ah 2008, 2010a, 2012a, 2013, 2014; Pan & Yuan 2012), and on the time-dependent 
hydrodynamical simulations of the model and their applications (e.g., Ruffert & Janka 1999; 
Lee et ah 2004; Setiawan et ah 2006; Lee et ah 2009; Carballido & Lee 2011; Sekiguchi & 
Shibata 2011; Caballero et ah 2012; Fernandez & Metzger 2013; Janiuk et ah 2013; Just et 
ah 2015). 

The structures and luminosities of the NDAFs in spherical coordinates had been studied 
(e.g., Liu et ah 2010a, 2012a, 2013). We stressed that the geometrical thickness of the disk 
was relevant to the advection and showed the existence of an empty funnel along the rotation 
axis which makes neutrino annihilation more efficient (Liu et ah 2010a). Furthermore, the 
distribution of the Bernoulli parameter, neutrino trapping, and components features, i.e., 
nucleosynthesis near the surface of the disk, were studied in the vertical direction (Liu et 
ah 2012a, 2013). The possible outflow might occur at the outer region of the disk, and the 
heavy nuclei tended to be produced near the surface of the disk, especially for the outer 
region. We also argued that the self-gravity effect may be important to the structure of the 
NDAF, and the neutrino luminosity further (Liu et ah 2014). The gravitational instability 
of the disk was reviewed, which might account for the X-ray flares or extended emission in 
GRBs. 

Convective motions have been widely studied in the black hole accretion disks. The verti¬ 
cally convective motions had been considered especially for the radiation pressure-dominated 
accretion disks (see, e.g., Bisnovatyi-Kogan & Blinnikov 1977; Shakura et ah 1978; Milsom 
et ah 1994; Blaes & Socrates 2001; S§.dowski et ah 2009, 2011; Gu 2012, 2015; Kawanaka & 
Kohri 2012). If the convective energy transfer in the vertical direction becomes effective, the 
inwardly advective energy would be suppressed. The convective motion in an accretion disk 
becomes important when the radiation pressure is dominant. As we know, the typical mass 
accretion rate of NDAF is about 0.1 — 1 M 0 s _1 , and the outer region of the NDAF should 
be radiation pressure-dominated if the convection is not taken into account (e.g., Liu et ah 
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2007; Kawanaka & Mineshige 2007). In the innermost region, the neutrino cooling becomes 
dominant, and the pressure would be dominated by gas. Thus we can study the convection 
motions in the NDAF directly. Kawanaka & Kohri (2012) first investigated the effects of 
the vertical convection in the NDAF model. If the vertical convection is considered, the 
hyperaccretion disk is expected to be hotter and the neutrino emission due to the electron- 
positron annihilation becomes the most efficient cooling process. They also presented that 
the sequence of the thermal equilibrium solutions for the convective NDAF had a viscously 
unstable branch, especially for the relatively small viscosity parameter, which might explain 
the origin of the drastic variabilities in the prompt emissions of GRBs. Recently, Jiang et 
al. (2014) investigated the super-Eddington accretion flows onto black holes using a global 
three-dimensional radiation magneto-hydrodynamical simulation. They found that the ra¬ 
diation of slim disk (e.g., Abramowicz et al. 1988) can be greatly increased via the vertical 
advection of radiation caused by magnetic buoyancy. As the natural extension of the slim 
disk in the state for very high accretion rate (Liu et al. 2008), NDAF should be expected to 
radiate a larger number of neutrinos due to the same mechanism. 

In this paper, we introduce the vertical convection in the thermodynamics of the NDAF 
model, and present its effects on the structure and luminosity of the hyperaccretion disk 
in spherical coordinates. The basic equations, boundary condition, and convective effect 
are presented in Section 2. In Section 3, the corresponding results of the structure, neu¬ 
trino luminosity and annihilation luminosity, and comparisons with the cases excluding the 
convection are shown. Summary is made in Section 4. 


2. Equations 
2.1. Basic equations 


As shown in Figure 1, a steady state axisymmetric accretion flow in spherical coordinates 
(R, 9, 0), i.e., d/dt = d/dcj) = 0 is considered (e.g., Kato et al. 2008). We adopt the 
Newtonian potential since it is convenient for self-similar assumptions. The basic equations 
of continuity and momentum are the following (see, e.g., Xue & Wang 2005; Gu et al. 2009): 
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where Vr, Vg, and v<p are the three components of the velocity, and p and p are the pressure and 
density of the disk. Here, the i?0-component of the viscous stress tensor are considered, i.e., 
Tr<p = pvRdivtj,/R)/dR. The kinematic coefficient of viscosity is written as u = ac 2 /f2 k, 
where the sound speed c s is defined as c 2 = p/p, the Keplerian angular velocity is fix = 
(GM/i? 3 ) 1 / 2 , M is the mass of the black hole, v$ = R£l k, and a is a constant viscosity 
parameter. Following Narayan & Yi (1995), the radial self-similar assumption is 
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According to substitute the above equations into Equation (l)-(4), which are trans¬ 
formed to the ordinary differential equations with respect of 9. Thus the continuity equation, 
i.e., Equation (1), is an identical equation, but we can define the mass accretion rate M to 
replace it as 

pit / 2 

M = —AnR 2 / pvRsm9d,9, (10) 

Je 0 

where 9o is the polar angle of the surface as shown in Figure 1, and the other physical 
quantities with subscripts “0” represent the values at the equatorial plane of the disk. 

Furthermore, the momentum equations, i.e., Equation (2)-(4), can be simplified as (e.g., 
Gu et al. 2009; Liu et al. 2010a, 2012a) 
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The equation of state is 


P 


Pgas + Prad + Pe + Pu 


pkT 1 + 3X nuc 
m p 4 


11 

+ I2 ar ‘ + 


2nhc 

3 


31 eP \ /G 
87rm p 3 


(14) 


where p gas , p ra d, p e , and p^ are the gas pressure from nucleons, the radiation pressure of 
photons, the degeneracy pressure of electrons, and the radiation pressure of neutrinos, re¬ 
spectively. k , T, m p , a, and h are the Boltzmann constant, temperature of the disk, proton 
mass, radiation constant and Planck constant, respectively. The equation can be applied to 
anywhere on the disk to describe the pressure, but we also need an equation to completely 
give the vertical distributions of the density, temperature, and pressure. For simplicity, we 
assume the polytropic relation, i.e., p = A'p 4 / 3 , in the vertical direction, K is a constant 
(e.g., Lee et al. 2005; Liu et al. 2008, 2010a, 2012a, 2013). Here we also assume the electron 
fraction Y c is 0.5 (e.g., Di Matteo et al. 2002; Gu et al. 2006; Liu et al. 2014), because its 
influences are very limited on the structure and neutrino luminosity of the disk, unless we 
care about the flavors of neutrinos emitted from the disk (e.g., Liu et al. 2013). X nuc is the 
mass fraction of free nucleons, and its expression obtained by taking Y e = 0.5 is (e.g., Kohri 
et al. 2005) 


-Xnuc = min[l, 295.5p 10 3 ' /4 T 1 9 / 8 exp(—0.8209/T n )], 


(15) 


where pio = p/10 10 g cm 3 and Tn = T/10 11 K. p u is the energy density of neutrinos, which 
can be expressed as (e.g., Lee et al. 2005; Liu et al. 2012a) 
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where T Ui = T a)Ui +T S:Ui is the total optical depth for neutrinos, T cltUi and are the absorption 
and scattering optical depth for neutrinos, and the subscript i runs for the three species of 
neutrinos z/ e , /y /t , and z/ T , which can be defined as (Di Matteo et al. 2002; Lee et al. 2005; 
Liu et al. 2012a) 
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where q Ui is the summation of the cooling rates per unit volume due to the neutrino reactions 
for different neutrinos, especially the Urea processes and electron-positron pair annihilation 
(e.g., Di Matteo et al. 2002; Liu et al. 2007). The neutrino scattering processes includes 
scattering by free protons, free neutrons, ct-particles, and electrons (e.g., Di Matteo et al. 
2002 ). 

The energy equation is written as 

Qv is Qadv T Qvi (19) 

where Q vis , Q a d v , and Q u are the viscous heating, advection and neutrino cooling rates per 
unit area, respectively (e.g., Liu et al. 2010a, 2012a, 2013). The heating due to the nuclear 
reaction can be neglected, because the gravitational energy is more powerful and efficient in 
the accretion systems, even for the hot dense matters in the NDAF. 

Q vis and Qadv can be expressed as 

r U 2 

Q vis = 2R q vis sin Odd, (20) 

Je o 


rn/2 

Q adv = 2R q adv sin Odd, (21) 

Je o 

where q v j S and q a d v are the viscous heating and advective cooling rates per unit volume. 
The viscous heating and advective cooling rate per unit volume are expressed as q v j S = 
vpR 2 [d(y < ] > /R)/dR] 2 and q a d v = pv R [de/dR— (p/p 2 )dp/dR\ (e is the internal energy per unit 
mass), respectively. After self-similar simplification, the expression is (Liu et al. 2010a, 
2012a) 


_ 9 api’l 
“ 4A 2 fV 


( 22 ) 


3 (p~ Pe)v R 



here the entropy of degenerate particles is neglected. 
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The neutrino radiation cooling rate Q v can be defined as (Lee et al. 2005; Liu et al. 
2012a, 2013) 
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As mentioned, different types of neutrinos are involved into different types of reactions and 
have their own optical depths. Additionally, we should point out that the prescription of the 
energy equation with convection is enough to describe the vertical energy transport (Lin et 
al. 2010a; Pan & Yuan 2012). 

Additionally, we follow the solutions of the boundary condition in Liu et al. (2012a, 
2013). As well as the principle of the Eddington luminosity, based on the balance between 
the vertical component of the gravity and radiation pressure at the surface of the disk, the 
surface temperature T s can be defined as 


= 


, 6GMm v 
acr T i? 2 


cot^n 


(25) 


where m u is the mean mass of a nucleon, ax is the Thompson scattering cross, and a is the 
radiation density constant. This equation is used to determine the angle of the disk surface 
and complete the basic equations, i.e., Equations (10)-(14) and (19). 


2.2. Vertical Convection 

In the vertical direction of the disk, the convective motion can carry the photons (or 
neutrinos) from the equatorial plane to the surface of the disk by the magnetic buoyancy. 
For slim disk model, the key issue is whether the vertical radiation transports process due 
to the magnetic buoyancy is faster than the radial advection process, thus the photons can 
escape from the surface before being advected into the black hole. In Jiang et al. (2014), the 
results indicated that the radiation can be widely increased because of the vertical advection 
of radiation caused by magnetic buoyancy. The average advective energy transport velocities 
along the vertical and radial directions at each radius had been calculated, which are with 
nonzero values and analogous to the bulk velocities, such as vr, c s or (e.g., S§.dowski et 
al. 2011; Jiang et al. 2014). As mentioned above, should NDAF be expected to radiate a 
larger numbers of neutrinos due to the same mechanism? 

Following the discussion in Kawanaka & Kohri (2012), we adopt the maximal value of 
the convective velocity physically allowed. Then we also estimate the typical timescale for 
the convective motion of the blob in the vertical direction t c and the advection timescale f a d v - 
The vertical convective speed along the vertical direction can be approximatively written as 

v c ~—vrcos9. (26) 

Here the components of the bulk velocity vr are cautiously chosen to present the radial 
advection and vertical convection (or advection), because there is no references to describe 



vertical convection well and we further consider the dynamical requirement. This is the 
vertically “infalling” speed onto the equatorial plane due to the vertical component of the 
gravity of the black hole. In other words, the convective velocity should be the same order 
of magnitude as the vertical components of vr, otherwise, the convection may be destroyed 
by the gravity and cannot exist in the vertical direction. Thus the convective timescale is 
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where R g = 2 GM/c 2 is the Schwarzschild radius. 


The vertical convection should exist if t c < t , adv , i.e., the polar angle of the convective 
region should satisfy 6 d < 6 < 9 C , where the critical polar angle 9 C is dehned by 
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When the vertical convection is included in the NDAF model and < 9 Cl Equation 
(21) can be replaced by 


rn/2 

Q adv = 2R q adv sin OdO, (30) 

Je c 

otherwise, the equation remains unchanged. Although only one equation is replaced, the 
essentials of the question has changed, i.e., how much ability of the vertical convection could 
suppress the radial advection and increase the vertical radiation transports? 


3. Numerical results 

The system of equations is closed, as Equations (10)-(14) and (19) with the boundary 
condition and given constant parameters R , M, M, and a. The main difference between 
the cases including and excluding the vertical convection is the description of the advection 
cooling rate. In our calculations we take M = 3M 0 , and define the dimensionless radius and 
mass accretion rate as r = R/R g and m = M/M & s -1 , respectively. 



9 

3.1. Structure 

Figure 2 (a) and (b) show the variations of the density and temperature at the equatorial 
plane with the dimensionless radius r, for which the given parameters are rri — 1 and a = 0.1. 
The solid and dashed lines correspond to the cases including and excluding the vertical 
convection, respectively. We notice that the variations of the density and temperature in 
the case excluding the vertical convection are similar to the solutions in Liu et ah (2010a, 
2012a, 2013). Considering the convection, the density slightly decreases and the temperature 
conversely increases, which are consistent with the results in Kawanaka & Kohri (2012). 
If the hyper accretion matter becomes hotter, the thickness of the disk is expected to be 
geometrically thicker due to the the vertical polytropic distribution of the temperature. 

Figure 3 shows the variations of the half-opening angle ( n/2 — 9 ) with the dimensionless 
radius r for the different cases. The black, red, and blue lines describe the cases that the 
given parameters (m, a) is (1, 0.1), (1, 0.05), and (0.1, 0.1), respectively. The solid and 
dashed lines correspond to (n/2 — 9 q) and (n/2 — 9 c ), and the convective motions exist in the 
area between these two lines to replace the radial advection. The thick and thin black solid 
lines correspond to the half-opening angle of the disk in the cases including and excluding 
the vertical convection. We also notice that the half-opening angle of the disk in the case 
excluding the vertical convection are similar to the solutions in Liu et al. (2012a, 2013), 
which has the slight difference comparing with the case including the convection. 

As the solutions in Liu et al. (2007, 2010a), the advection and neutrino cooling mech¬ 
anisms are respectively dominant cooling modes in the outer and inner region of the disk. 
Thus the effects of the vertical convection mainly reflect in the inner region, i.e., 9 C is closed 
to 7 t/ 2, until very near the black hole, R < 10 R g . The amplitudes of the differences of the 
density and temperature between the cases including and excluding the convection can also 
be understood by the above explanations. 

Besides, we also notice that the accretion rate and viscous parameter can markedly 
change the polar angle 9 0 and 9 C , and the effects of the viscous parameter are much more 
significant. Higher density and temperature is conducive to the more neutrinos emission, and 
the density and temperature decrease gradually from the equatorial plane to the surface of 
the disk, so the number of neutrinos radiated from the equatorial plane of the disk are much 
more than that from other regions, especially nearby the surface. The convective energy 
transfer in the vertical direction can be effective to suppress the radial advection in the 
NDAF, but the region dominated by the vertical convection is deviated from the equatorial 
plane of the disk, thus it can be expected that the neutrino emission rate should increase 
slightly. 
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3.2. Neutrino Luminosity and annihilation luminosity 


Using the neutrino cooling rate Q„, the neutrino luminosity before annihilation L u is 
obtained as 


u 


r»300 R a 


47T 


Q u RdR. 


13 R a 


(31) 


Furthermore, in Liu et al. (2010a), the luminosity of neutrino annihilation L vD can be 
roughly evaluated by the assumption: rj oc V~^ (see, e.g., Mochkovitch et al. 1993), where 
rj = L V v/L v is the annihilation efficiency, and V ann is the volume above the disk. If we take 
the outer boundary of the disk R ont is 300i? 3 , U ann can be estimated by integrating the region 
of 9 < 9o and R < R 0ut . We can estimate the proportionality constant of the efficiency to the 
volume above the disk as a function of accretion rate from the vertically-integrated NDAF 
model where 9q corresponds to n/2 (e.g., Liu et al. 2007). 


Figure 4 displays the neutrino luminosity L y and annihilation luminosity L uD for the 
varying dimensionless mass accretion rates m. The solid and dashed lines correspond to the 
cases including and excluding the vertical convection, respectively. The thick and thin lines 
display the estimations of the neutrino luminosity and annihilation luminosity, respectively. 
The luminosity has a certain increase, which is origin from the suppressed advection and 
correspondingly increased neutrino cooling in all radii of the disk. The influence on the inner 
region of the disk is stronger in the cases of the lower accretion rate (e.g., m = 0.1) than in 
the cases of larger accretion rate (e.g., m = 1) as shown in Figure 3. Because neutrinos are 
mainly launched from the inner region of the disk, the difference of the luminosity between 
with and without the convection in the case of the lower accretion rate is slightly larger 
than that in the case of the larger accretion rate. Additionally, in Kawanaka & Kohri 
(2012), they found that the electron-positron pair annihilation replaces the Urea process as 
a dominant neutrino emission process when convection works efficiently. These two processes 
are intensely related to the density and temperature (e.g., Di Matteo et al. 2002; Liu et al. 
2007). In our calculations, because the density and temperature do not change apparently, 
the annihilation becomes stronger but still weaker than the Urea processes until for the larger 
radius. 


Considering with the increase of the opening angle, the neutrino annihilation luminosity 
will also be magnified as shown in the Figure 4. Obviously, the neutrino luminosity and 
annihilation luminosity derived from the calculations considering the effects of the opening 
angle of the disk are significantly larger than those derived from the calculations neglecting 
the thickness of a disk as shown in the above estimation method of annihilation luminosity 
(e.g., Liu et al. 2010). We notice that for M ~ 5 M 0 s -1 , the density of radiated neutrino 
is so large that the annihilation efficiency is close to 1. Furthermore, the disk is extremely 
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thick, which restrains the ejection existing in a vary narrow empty funnel along the rotation 
axis. As discussed in Liu et al. (2010a), since the large opening angle of the disk and 
the sufficient output of energy, the neutrino annihilable ejection required by GRBs may be 
naturally explained. 

4. Summary 

For one of the radiation pressure-dominated accretion disks, we study the effects of the 
vertical convection on the structure and luminosity of the NDAF around a stellar-mass black 
hole in spherical coordinates. The major points we wish to stress are as follows: 

1. The convective energy transfer in the vertical direction can be effective to suppress the 
radial advection for the different accretion rates and viscous parameters in all radii of 
the NDAF, but the region dominated by the vertical convection is deviated from the 
equatorial plane of the disk. 

2. The density, temperature, and opening angle of the disk are slightly different between 
two cases with and without the vertical convection. 

3. Considering the vertical convection, the neutrino luminosity and annihilation luminos¬ 
ity have a certain increase. The increased opening angle of the disk can enhance the 
annihilable efficiency and the collimation of the ejecta, which further verifies that the 
NDAF can be regarded as the candidate of the central engine of GRBs. 

In future, we should work on the numerical calculations or simulations of the NDAF 
model with more precise description of the vertical convection and the neutrino transfer to 
replace the rough estimates by the convective region and polytropic relation in this work. 
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the CAS Open Research Program of Key Laboratory for the Structure and Evolution of 
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Fig. 1.— Schematic picture of the system composed of a black hole and an accretion disk in 
spherical coordinates, similarly see also in Liu & Xue (2011). 
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Fig. 2.— Variations of the density and temperature at the equatorial plane (p 0 and T 0 ) with 
the dimensionless radius r, for which the given parameters are m = 1 and a = 0.1. The solid 
and dashed lines correspond to the cases including and excluding the vertical convection, 
respectively. 
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Fig. 3.— Variations of the half-opening angle (vr/2 — 9) with the dimensionless radius r for 
the different cases. The solid and dashed lines correspond to (vr/2 — 9 0 ) and (vr/2 — 9 C ), 
respectively. The thick and thin black solid lines correspond to the half-opening angle of the 
disk in the cases including and excluding the convection with m = 1 and a = 0.1. 
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Fig. 4.— Neutrino luminosity L v and annihilation luminosity L vU for varying dimensionless 
mass accretion rates fn. The solid and dashed lines correspond to the cases including and ex¬ 
cluding the vertical convection, respectively. The thick and thin lines display the estimations 
of the neutrino luminosity and annihilation luminosity, respectively. 




